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A microscopic model for the room-temperature multiferroic BiFeOa that includes two 
Dzyaloshinskii-Moriya interactions and single-ion anisotropy along the ferroelectric polarization pre¬ 
dicts both the zero-field spectroscopic modes as well as their splitting and evolution in a magnetic 
field. Due to simultaneously broken time-reversal and spatial-inversion symmetries, the absorp¬ 
tion of light changes as the magnetic field or the direction of light propagation is reversed. We 
discuss three physical mechanisms that may contribute to this absorption asymmetry known as 
non-reciprocal directional dichroism: the spin current, magnetostriction, and single-ion anisotropy. 

We conclude that the non-reciprocal directional dichroism in BiFeOs is dominated by the spin- 
current polarization and is insensitive to the magnetostriction and easy-axis anisotropy. With three 
independent spin-current parameters, our model accurately describes the non-reciprocal directional 
dichroism observed for magnetic field along [1,—1,0]. Since some modes are almost transparent 
to light traveling in one direction but opaque for light traveling in the opposite direction, BiFeOs 
can be used as a room-temperature optical diode at certain frequencies in the GHz to THz range. 

Our work demonstrates that an analysis of the non-reciprocal directional dichroism spectra based 
on an effective spin model supplemented by first-principles calculations can produce a quantitative 
microscopic theory of the magnetoelectric couplings in multiferroic materials. 

PACS numbers: 75.25.-j, 75.30.Ds, 75.50.Ee, 78.30.-j 


I. INTRODUCTION 

BiFeOs is the only material known to exhibit multi¬ 
ferroic behavior at room temperature. Because its fer¬ 
roelectric (FE) transition temperaturei Tc ~ 1100 K is 
significantly higher than its Neel transition temperature^ 
Tn ~ 640 K, BiFeOs is a type I multiferroic. Although 
the non-magnetic FE polarization^ ~ 90 /iC/cm^ is 
much larger than the magnetic contribution^-— induced 
by the distorted spin cycloid^i^-— , the magnetic domain 
distribution of BiEeOs can be manipulated by an applied 
electric field^ii^. 

A great deal has been learned about BiEeOs since 
the first single crystals became available for inelas¬ 
tic neutron scatteringi^"— , Raman scatteringi^iii, and 
THz spectroscop} d^d^ measurements. It is now under¬ 
stood that two sets of interactions control the cycloid 
of BiEeOs: two Dzyaloshinskii-Moriya (DM) interactions 
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produced by broken inversion symmetry and a single-ion 
anisotropy2£ (ANI) K along the direction of the EE po¬ 
larization Whereas the DM interaction^ Di per¬ 

pendicular to is responsible for the long 62 nm cy¬ 
cloidal period, the DM interaction^i^-— D 2 along P^^ 
is responsible for a small cycloidal tilt^. Above the crit¬ 
ical magnetic field the cycloidal tilt develops into 
the weak ferromagnetic (EM) moment^i^i^i of a G-type 
antiferromagnet (AE) that is isosymmetrically canted 
by an antiferrodistortive rotation (R^[1,1,1]) of the R3c 
structure^^. 

Inelastic neutron scattering measurementsi^-— were 
used to extract the AE nearest- and next-nearest neigh¬ 
bor exchange interactions^^ Ji = —5.32 meV and J 2 = 
—0.24 meV between the S = 5/2 Ee^+ spins on the 
pseudo-cubic unit cell sketched in Eig.l(a) with lat¬ 
tice constant a = 3.96A. However, those measure¬ 
ments lacked the sensitivity to resolve the ordering 
wavevectors on either side of the G-type AE wavevec- 
tor Qo = (27r/a)[0.5, 0.5,0.5] at (27r/a)[0.5 ± (5, 0.5, 0.5 ^ 
(5], where S ~ 0.0045 is inversely proportional to the 
cycloidal period a/{\/25). Recent neutron scattering 
measurements^ with higher precision were able to distin¬ 
guish the two cycloidal ordering wavevectors and found 
thati2^ Di = 0.18 meV and K = 0.0039 meV. But even 
those measurements lacked the precision to obtain D 2 , 
which was set to zero. 

By contrast, the frequencies of the spin-wave (SW) 
modes at the ordering wavevector Q can be pre¬ 
cisely measured with Raman scatteringi^ and THz 
spectroscop} d^d^ . The parameters K, Di, and D 2 were 
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estimated by fitting the frequencies^i of the four observed 
zero-field THz modes. With no remaining adjustable pa¬ 
rameters, that same model predicted^ the evolution and 
activation of the THz modes^^ in a magnetic field along 
[ 0 , 0 , 1 ]. 

We now use this microscopic model to predict the 
asymmetry A(a(co’) in the absorption a{uj) of light when 
the direction of the magnetic field or, equivalently, 
the direction of light propagation is reversed. Called 
non-reciprocal directional dichroism (NDD), absorption 
asymmetry was first observed by Hopfield and Thomas^ 
over 50 years ago in CdS. Much more recently, the pre¬ 
cise symmetry requirements for NDD in magnetic ma¬ 
terials were systematically investigated by Szaller et al. 
[s^. Strong NDD is expected for the spin excitations of 
multiferroic materials when both time reversal and spa¬ 
tial inversion symmetries are broken by the spin state. 
Both the magnetic and electric components of THz ra¬ 
diation can excite SWs in multiferroic materials. The 
NDD exhibited by simultaneously electric- and magnetic- 
dipole active excitations has been extensively studied in 
Ba2CoGe207 [111-0, Sr2CoSi207 H, Ca2CoSi207 H, 
Gd-o.sTbo.sMnOa [ssj . and Euo. 55 Yo. 45 Mn 03 [36| . 

Because the cycloidal spin state is produced by the 
competition between DM, exchange, and ANI interac¬ 
tions, three distinct physical mechanisms can produce 
NDD in BiFeOs: the spin current (SC) driven by the 
DM interactions, magnetostriction (MS) or the electric- 
field induced changes in the exchange interactions, and 
the electric-field induced changes in the ANI. Remark¬ 
ably, the dynamical magnetoelectric coupling governing 
the NDD in BiFeOs is dominated by the two sets of 
SC polarizations associated with Di and D 2 . Qualita¬ 
tively, the SC dominates the magnetoelectric coupling 
in BiFeOs because spin fluctuations (58^ are transverse 
to the almost collinear, cycloidal spin state ( 8 ^). Since 
^8i X (8j) 7 ^ 0 but (^8^ • {Sj) ^ 0 (for nearby sites i and 
j) and the ANI is extremely weak, spin fluctuations more 
strongly affect the SC-induced polarization than the MS- 
and ANI-induced polarizations. 

As a fraction of the total light absorption at a given 
wavelength, NDD is most pronounced for a mode with 
fluctuations out of the cycloidal plane at 15.5 cm“^. At 
this wavenumber, BiFeOs is almost transparent for light 
traveling in one direction but opaque for light traveling in 
the opposite direction. Therefore, BiFeOs can be used as 
an optical diode that operates up to room temperature. 

Despite the success of our model describing the NDD 
for magnetic field along [ 1 ,— 1 , 0 ], several questions re¬ 
main open. Although our model predicts NDD to be 
absent for light propagating along k = [ 0 , 0 , 1 ], a static 
magnetic field along [ 7 ^, 7 ^,/ 1 :], and THz electric-field ori¬ 
entation e = [1,1,0] or [1,—1,0], weak NDD has been 
observed for a magnetic field along [ 1 , 1 , 0 ] under those 
conditions. An optical misalignment of the THz electric- 
and magnetic-field vectors e and h = k x e may be re¬ 
sponsible for this effect. In addition, the mean absorption 
a{uj) (the absorption a{uj) averaged over positive and 




FIG. 1: (Color online) (a) The exchange interactions Ji and 
J 2 on the pseudo-cubic lattice for BiFeOs with Fe^^ ions at 
the corners of the cube and two hexagonal layers shown. Due 
to the rhombohedral distortion along z', J 2 and J 2 may be 
slightly different, (b) The spin state in zero magnetic field, 
both with electric polarization along z^ The canting of the 
spins out of the {x plane is indicated by the angle r in 
the inset. The variation of the canted component along x' is 
also shown by the dashed line in (b). 

negative magnetic fields) is not as accurately predicted 
by our model as the NDD. 

This paper complements a recent worb^i that presents 
detailed experimental results for both the individual ab¬ 
sorptions and the NDD. We have divided this paper into 
six sections. Section H presents our microscopic model 
and Section HI presents the predicted mode frequencies. 
Section IV describes the three polarization mechanisms 
and presents results for the magnetization and polariza¬ 
tion matrix elements, with symmetry relations provided 
by Local Spin-Density Approximation (LSDA)-^!/ cal¬ 
culations. Results for the NDD are presented in Sec¬ 
tion V. Section VI contains a discussion and conclusion. 
While Appendix A summarizes the experimental details. 
Appendices B, C, and D treat the SC-, MS-, and ANI- 
induced polarizations, respectively. For convenience, the 
unit vectors used in this paper are given in Table I. 

II. MICROSCOPIC MODEL 

In a magnetic field H = iLm, the spin state and SW 
excitations of BiFeOs are evaluated from the microscopic 
Hamiltonian 

^ = - Ji ^ s, • s,- - J2 ^ Si • s,- 

+Di y](z' X eij/a) • (S* x Sj) 

+D 2 ^ (-1)”^ z' • (S, X S,) 

(di) 

i i 

where eij = ax, ay, or az connects with its nearest 
neighbor Hj = R^ + Bij. Since the unit vector z' points 
along a cubic diagonal parallel to the FE polarization 
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TABLE I: Unit vectors 


X, y, Z 

Pseudo-cubic laboratory reference frame 

Z = Z^ 

Orientation of the electric polarization 
along one of the cubic diagonals 

x', y', z' 

Cycloidal reference frame 

x' v' 7' 

1 J mi ^m 

Cycloidal reference frame for domain m 

U 

X, y, or z 

m 

Orientation of the static magnetic field 

Hi 

Local single-ion ANI axis 

k 

Direction of light propagation 

e 

Orientation of the THz electric field 

h 

Orientation of the THz magnetic field 


the Di sum has the form proposed by Katsura et 
alM. The hexagonal layers normal to z' are separated 
by c = a/v^ and are labeled by the integer • 

z'/c. Consequently, the D 2 sum alternates sign from one 
hexagonal layer to the next. Notice that the local DM 
interactions Di (z! x Bi^jja) and D 2 z' are, respectively, 
perpendicular and parallel to zl. 

There are eight possible orientations for || z! 

along the four cubic diagonals. For every possible z', the 
three magnetic domains have different x' and y'. When 
zl = [1,1,1] (all unit vectors in Table I are assumed nor¬ 
malized to 1 ), the possible orientations for the x' axis 
are x'l = [ 1 ,- 1 , 0 ], X 2 = [ 1 , 0 ,- 1 ], and Xg = [ 0 , 1 ,- 1 ] 
with corresponding = z' x These three magnetic 
domains have cycloidal ordering wavevectors 


Qm 


_ 2^/2tt5 , 

Qo H- 


( 2 ) 


Hence, the ordering wavevectors for different domains are 
Qi = (27r/a) [0.5 + (5,0.5 - (5,0.5], Q 2 = (27r/a)[0.5 + 
(5,0.5,0.5 - 5 ], and Q 3 = (27r/a)[0.5,0.5 + (5,0.5 -5]. In 
terms of ^ ^ 1 , the period of the cycloid in zero field is 
a/ (a/ 2(5) ^ 62 nm. 

As mentioned above, the DM interactions Di and 
D 2 only couple nearest-neighbor sites. In a previous 
formulatioi> 2 ii 2 S of this microscopic model, Di coupled 
next-neighbor sites within the same hexagonal layer. Due 
to the very long cycloidal period p ^ a of BiFeOs, 
the equilibrium and dynamical properties of these two 
Hamiltonians are the same up to errors of order ~ 
2 X 10“^. Specifically, earlier predictions for the SW 
mode frequencies^i^ and critical magnetic field^ are 
unchanged. However, the earlier DM interaction Di is 
now multiplied by ^/2. Because the nearest-neighbor DM 
interactions are much larger than those between next- 
neighbor spins, the Hamiltonian above provides a close 
connection with recent first-principles calcnlation a^^d^ . 

Since the Di and D 2 terms in H depend only on 
z', H is independent of the magnetic domain. For a 
specific domain m, the first SC term can be written 
ysc _ y^DiNy' • Ti, where 


Ti = ^ y] |Si X Sjj, TO = 1 , 2 , (3) 

(ij)- 


= ^ E{s*xSi}, m = 3, (4) 

ihj)^ 

where is a sum over nearest neighbors with Hj — 

Hi = au. These relations assume that the spins on 
each hexagonal layer depend only on the integer r = 
a/2x' 'Hi/a. So for domain 2, S(Ri + ax) = S(Ri — az). 
The cross products in Eqs.(|3]) and (|1]) couple spins with 
indices r and r -1- 1 on neighboring layers. 

The second SC term ¥ 2 ^ proportional to D 2 can be 
written 1 ^ 2 ^^ = D 2 N z' • T 2 , where 

T2 = ^E(-1)”’{s*xS,}. (5) 

(*d> 

Like ¥ 2 ^ also couples neighboring spins on neigh¬ 

boring layers. 

The nearest- and next-nearest neighbor exchange 
interactions^^ Ji = —5.32 meV and J 2 = —0.24 meV 
were obtained from recent inelastic neutron scattering 
measurement si^“— between 5.5 meV and 72 meV. On 
the other hand, the small interactions Di, D 2 , and K 
that control the cycloid can be obtained from THz spec¬ 
troscopy measurement si^i21 below 5.5 meV (44.3 cm“^) 
in zero magnetic field. 

We have neglected the broken spatial symmetry be¬ 
tween the exchange interactions due to the rhombohe- 
dral distortion. While all Ji interactions must remain 
the same due to the rotational C 3 symmetry about z', J 2 
may reflect the rhombohedral distortion. For example, 
next-nearest neighbors separated by Ri = a(x -f y) and 
R 2 = a(x — y) may experience slightly different exchange 
interactions, denoted by J 2 and J 2 in Fig. 1(a), because 
Ri • z' = 2a/v /3 while R 2 • z' =0. However, based on the 
excellent agreement between theory and experiment for 
the mode frequencies reported in Section HI and because 
J 2 is already so small, we expect this exchange anisotropy 
to have a very minor effect on the NDD. 

For a given set of interaction parameters, the spin state 
of BiFe 03 is obtained by minimizing the energy E = {H) 
over a set of variational parameters'^. Fixing 6 = 1/q^ 
where g ^ 1 is an integer, the energy E is minimized over 
the variational parameters on a unit cell with q sites along 
x' and two hexagonal layers. The spin state on layer n is 
assumed to be identical to the spin state on layer n -h 2. 
The wavevector parameter 6 is determined as a function 
of field by an additional minimization loop over q. In 
zero field, 5 ^ 0.0045 and q = 222. We verify that the 
corresponding spin state provides at least a metastable 
minimum of the energy E by checking that the classical 
forces on each spin vanish. 

Ignoring the cycloidal harmonics C/>i produced by D 2 
and K but including the tilt^ r produced by D 2 , the spin 
state in zero field can be approximated by 



= 5'(—I)’^^^ COSTsin(27r(5r). 

( 6 ) 

SyiK) 

= 5'sinr sin(27r(5r). 

( 7 ) 

S,'(R) 

= 5'(—I)’^^^ cos(27r(5r). 

( 8 ) 
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FIG. 2: (Color online) The predicted critical field (solid) ver¬ 
sus So for field orientation m = [0, 0,1]. The horizontal line is 
the experimental valued for He- Also plotted is the single-ion 
ANI K (dash-dot) versus So- 

This tilted cycloid is plotted in Fig. 1(b). Averages 
over this state are readily performed using {Six'‘^) = 
(5'^/2) cos^ r, {Siy'‘^) = (5'^/2) sin^ r, and {Siz'‘^) = 
S‘^/2. In zero field, averages over the tilted cycloid are 
fairly accurate because^i even harmonics like C 2 vanish 
and Cs ~ 5 x 10“^. Corrections to the averages are then 
of order « 2.5 x 10“^. 

For comparison, the spin state of the canted AF at zero 
field can be simply written in terms of the canting angle 
r within the {x'^y'^z'} coordinate system as 

Sn = 5'[(— COST, sinr, 0] (9) 

on hexagonal layer n. Recall thafi^i sinr = So/S where 
2iibSo is the weak FM moment of the AF phase along 
y' above He- Whereas susceptibility measurements^ indi¬ 
cate that 5'o = 0.015, a recent neutron-scattering study^i 
suggests that So ^ 0.05 equivalent to r ^ 1°. By con¬ 
trast, LSDA+1/ (1/ = 5 eV)^ gives So = 0.014, in agree¬ 
ment with the former experimental result. Note that 
D 2 = — 2Ji5'o/*S' = — 2Jisinr is a linear function of So 
and of sinr ~ r. 

We now adopt a different approach to estimate D 2 . 
The three parameters Di, 1^2, and K are fixed by two 
conditions: the period of the cycloid must match the 
measured period and and the frequencies of the four pre¬ 
dicted SW modes in zero field must match the measured 
frequencies^^. A third condition is provided by the de¬ 
pendence of the predicted critical field He on So- As 
shown in Fig.2, the measured critical field of 18.8 T for 
m = [0, 0,1] requires that2^ So = 0.02, corresponding to 
r = 0.008 or 0.45°. While Di ^ 0.180 meV is virtually 
independent of 5'o, D 2 linearly increases with So- Figure 
2 indicates that K increases almost quadratically with 
So from a value of K = 0.0031 meV when 5'o = 0. Cor¬ 
responding to So = 0.02, we obtain D 2 = 0.085 meV and 
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FIG. 3: (Golor online) The mode spectra at multiples of the 
ordering wavevector in an extended zone scheme (a) without 
higher harmonics of the spin and (b) with higher harmonics 
for the predicted parameters of BiFeOa- 


K = 0.0051 meV. A somewhat smaller value K = 0.0039 
meV was given in Ref.[26[, which took 1^2 = 0. 

With other parameters fixed and m = [1,—1,0], a 
value of D 2 smaller than about 0.079 meV would stabilize 
a different canted AF phase above He with spins tilted 
above and below the {x 'plane due to the dominant 
single-ion ANI. Hence, the coplanar AF phase of Eq.(|9]) 
is barely stabilized by the second DM interaction. 


III. SPECTROSCOPIC MODE FREQUENCIES 

Using the parameters given above, H predicts the evo¬ 
lution of the modes with magnetic fiel d^^i^^i^^ for all 
orientations m. SW modes at the ordering wavevec¬ 
tor Q can be labeled^ as in-cycloidal-plane modes 
and out-of-cycloidal-plane 4^^ modes. In an extended 
zone scheme, those mode frequencies are plotted ver¬ 
sus q/5 for wavevector (27r/a)[0.5 + 0.5, 0.5 — q] in 
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Fig.3(a). For simplicity, and denote both the 
modes and their frequencies. Neglecting higher spin har¬ 
monics, ^rn = 1^1 ^0 and It folloWS 

that <Fi = ^ 0 - 

Higher harmonics generated by the tilt and ANI split 

fl 2) (1 2) 

each mode with m > 1 into two labeled ^ or . 

For the predicted parameters of BiFeOs, those modes are 
plotted versus wavevector in Fig.3(b). While the m = 1 
modes are strongly affected by the spin harmonics, the 
former mode scheme remains fairly accurate for m > 1. 
Because the splitting of the low-frequency modes was not 
considered, recent Raman stndie a^^i^^ misidentified the 
observed modes with some out-of-plane modes mistaken 
for in-plane modes and vice versa. 

Despite the substantial splitting of and 
is only slightly larger than ^o- The nearly degener¬ 
ate and ^0 modes cannot be separated by THz 

measurementin zero field. 

In Fig.4, the predicted and measured^ mode frequen¬ 
cies are plotted versus field for orientations m = [0, 0,1], 
[1,1,0], and [1, —1, 0]. Experimental data was not avail¬ 
able for the THz modes above He for the last two field 
orientations. The experimental results for m = [1,1,0] 
and [1, —1, 0] are presented here for the first time with ex¬ 
perimental details summarized in Appendix A. Because 
its frequency was too low, ^ was not detected when 
m = [1,1,0] and [1, —1, 0]. The predicted mode frequen¬ 
cies of the stable domain(s) are presented in the solid 
curves: domain 1 for m = [0,0,1] and [1,1,0] and do¬ 
mains 2 and 3 for m = [1, —1,0]. For m = [0,0,1], the 
mode that dips below ^ arises from metastable do¬ 
mains 2 and 3, as seen by the agreement with the dashed 
curve. Hence, metastable domains may survive up to 
about 10 T. 

With So = 0.02, the agreement between experiment 
and theory is even better than previously reportedi^ for 
m = [0,0,1] with So = 0.015. Nevertheless, that agree¬ 
ment deteriorates somewhat above 12 T, particularly for 
m = [1, — 1, 0], when avoided mode crossings strongly af¬ 
fect the mode frequencies. It is possible that the trial spin 
state is not sophisticated enough at high magnetic fields. 
For example, the spin state in high magnetic fields may 
have a periodicity greater than two hexagonal layers. 

Above i^c, the canted AF state of Eq.(|9|) supports only 
two modes that are labeled a and [3 in Eig.4. Because the 
transition at He is first order, the spectroscopic modes 
change discontinuously at the critical field. 

The estimates given above for Di, and D 2 were 
based on fits to the four THz modes observe d^^d^ in zero 
field^. Experimental data points in Eig.4 indicate that 
those four modes correspond to (nearly degen¬ 
erate), and with frequencies 16.2, 20.7, 

22.4, and 27.6 cm“^, respectively. 



H(T) 


FIG. 4: (Color online) The theoretical mode spectra (solid 
curves) and experimental measurements (boxes) versus field 
for field orientations m = (a) [0,0,1], (b) [1,1,0], and (c) 
[1, — 1, 0]. Solid vertical lines make the transition to the canted 
AF state. The dashed curve in (a) indicates the predicted 
for metastable domains 2 and 3. 
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IV. POLARIZATION MATRIX ELEMENTS interaction: 


At zero magnetic field, only a few of the SW modes 
are optically active with finite magnetic-dipole resonance 
matrix elements (n|M|0), where 

( 10 ) 


pSC^-^^Wjx'xTi}^, (12) 

with yv = 0, and = 

-A()t/(2\/2) so that Ay = -Ay/2. 

The second SC-induced polarization is associated with 
the DM interaction D 2 '. 


is the magnetization operator, |0) is the ground state with 
no SWs, and \n) is the nth excited state with a single SW 
mode at the cycloidal wavevector Q. At finite magnetic 
fields, all SW modes also have non-zero matrix elements 
^^|pind|o) Qf induced electric polarization The 

coexistence of the magnetic-dipole and polarization ma¬ 
trix elements is responsible for the NDD observed in the 
THz absorption spectra for field along [1,—1,0]. The 
physical mechanisms that contribute to P^^^ below Tn 
can be divided into three classes: SC, MS, and ANL 
For the SC- and MS-induced polarizations, we use 
LSDA+I/ calculations^ to simplify the matrices con¬ 
necting the induced polarizations with the spin operators. 
This greatly reduces the number of polarization param¬ 
eters. In some instances, those matrices were simplified 
even further, either because some matrix elements were 
roughly equal or because additional matrix elements had 
a negligible effect on the NDD. Those additional simpli¬ 
fications are described in Ref. [io| . This section expresses 
the induced polarizations in the cycloidal reference frame 
{x'^y'^z'}. In the laboratory reference frame {x^y^z}^ 
the induced polarizations are given in Appendices B, C, 
and D. 


A. SC-induced polarizations 

The SC-induced polarization P^^ is produced by shifts 
in the O locations due to the hopping of electrons between 
Fe 3d and O 2p orbital s^^d^d^ . The first SC-induced 
polarization is created by the well-known inverse DM 
interaction^i^i^ corresponding to the Di term in the 
Hamiltonian. This polarization can be generally written 
as 


pSC 

^2a 




2a : 


(13) 


where T 2 was defined by Eq.(|5]). As shown in Appendix 
(2) 

B, the z' coefficient may differ from the x' and y' 
coefficients = A^,^^. 

X y 

For the simple tilted cycloid of Eqs. (03) in zero mag¬ 
netic field. 


(P^^) = 27r5'^(5cosr 





A^^) z' 


(14) 


(pSC) = sm2Tz'. (15) 

If the cycloid were not tilted, only the first SC polariza¬ 
tion would be nonzero. When c = d, = Ai!L/2 and 

z y y z I 

the first term in (Pi^) along y' vanishes. 


B. MS-induced polarizations 

The first MS-induced polarization is produced by the 
uniform displacement of Ee with respect to O: 



(16) 

= Cl,, z' ■ Wi, 

(17) 


(18) 


i'iJ) 


(11) 

/3 

where Ti was defined by Eqs. ([3]) and (|1]). According 
to Eg. fjBIip . the four nonzero matrix elements of A^^^ 
are A^^, = -A^, = ±(c - d), A^, = -2V2c, and 

= — \/2d, where the plus sign is for domain 2 and 
the minus sign is for domains I or 3. 

In a simplified version of the first SC-induced polariza¬ 
tion with c = d, the diagonal terms A^!^, and A^!^, would 
vanish. Then = X^^^Tiy^ and Pfjy = 

would reduce to the usual forn>2^ for the inverse DM 


It is easy to show that x' • Wi =0. Eor a simple twisted 
cycloid, 

^pMS^ _ —x/^Ciz'S‘^cos^rz'. (19) 

The energy —E • P^^ uniformly shifts all the nearest- 
neighbor interactions by AJi = C\z'Ez' j\f3. 

The second MS-induced polarization can be written^ 

P^^ = C 2 z'x W 2 , (20) 


( 21 ) 
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Unlike Wi, W 2 alternates sign from one hexagonal layer 
to the next. The cross product with z' in Eq. (| 2Q|) ensures 
that remains a polar vector—. For a simple tilted 
cycloid in zero field, (W 2 ) = (P^^) = 0 . The energy 
—E-P^^ shifts the nearest-neighbor exchange interaction 
Ji by an amount proportional to C 2 E. For example, 
the nearest-neighbor exchange between spins at and 

Ri + az is shifted by AJi = — Ey)/\/3. 

Appendix C shows that C 2 = V^Ciy'. 

The MS-induced polarization associated with next- 
nearest neighbor sites can be similarly constructed start¬ 
ing with 

= ^ E Si-S,-, (22) 

= ^ E (23) 

where all next-nearest neighbor pairs are double 

counted with Hj — R^ = av, |v| = \/2, and v • u = 0. 
So for u = X, V = (0,1,1), (0,1,—1), (0, — 1,1), and 
(0, — 1, — 1). For next nearest neighbors, both R^ and Rj 
lie on either even or odd layers. 


Since x'-Ws = x'-W 4 = 0 , the polarizations associ¬ 
ated with W 3 and W 4 are 

^t'E=C^3.4yy'-W3.4, 

(24) 

EA' = <^3,4.' Z' • W3,4. 

(25) 

For a simple twisted cycloid. 


(P“S) = 2^/lS^Czz' z' 

(26) 

while (P^®) = 0. The energy —E ■ Pg^® 

uniformly 


shifts all the next nearest-neighbor interactions by A J 2 = 

Csz'Ez' IV^. 

Another possible MS-induced polarization is associ¬ 
ated with the spin exchange ANI or different exchange 
couplings for different spin components Sia- Because it 
is of order this polarization can be neglected. 

C. ANI-induced polarizations 

The ANI-induced polarization P^^^ = pANi ^ pA^^, 
which arises from the spin-dependent hybridization be¬ 
tween the Fe ions and their ligands, contains compo¬ 
nents perpendicular or parallel to z'. As shown in 
Appendix D, the perpendicular polarization pANi _ 
pANi(i) ^ pANi(2) terms associated with 

the electric-field dependence of the local single-ion ANI 
axis [sin Oi cos , sin Oi sin , cos Oi] defined by Fq. m- 
The first set is produced by the dependence of the polar 
angle Oi on the electric field E: 

^ + Siy'y')Si.^, (27) 


which agrees with the first ANI-induced polarization pro¬ 
posed by deSousa et al. ji^ . 

An additional perpendicular polarization 

E E{(5»'' - Siy,^)y + (28) 

i 

with ^2 = Ci/(2a/ 2) was proposed in Ref.[48|. However, 
the cross terms Sij^Si^ {/3 7 ^ 7) in Eq. (jD5p cancel this 
contribution. 

The second set of perpendicular ANI-induced terms is 
produced by the dependence of the azimuthal angle (pi 
on E: 

pANi(2) ^ ^(-1)-^ (29) 

i 

which was not previously proposed. 

We also construct the ANI-induced polarization paral¬ 
lel to z' produced by the electric-field dependence of the 
constant K: 

i 

which shifts the single-ion ANI by AK = ^4Ez'/4:. For 
a simple tilted cycloid in zero field, (pANi^ _ /S 

includes only a contribution from ^4 and is parallel to z '. 


D. Total induced polarization 

With all proposed terms, the net induced polarization 
in the cycloidal phase is + P^^ + pANi^ 

the simple tilted cycloid, 

(z' ■ = S^\-2ttX[]15cost - E 7 sin2T 

+ - 73(7i^,cos2r + 2V3C3^.}. (31) 

Of course, the components of perpendicular to z' 

do not change the magnitude of the total polarization 
significantly. The change in polarization from the para¬ 
magnetic phase above Tn to the cycloidal phase below 
Tn is given by Eq. m- Recently, Lee et al [3] observed 
that has a magnitude of about 400 nC/cm^ and 

opposes P^^ due to the suppressed displacement of the 
Fe ions compared to the Bi ions. 

By comparison, the induced polarization of the canted 
AF evaluated using Eq.(|9]) is given by 

(z' • Paf) = —sin 2 r VsCiz' cos 2 r 

+2^Czz’], (32) 

which has no ANI contribution because the spins are in 
the {xP y'} plane. So the change in polarization from the 
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AF phase to the cycloidal phase at zero field is given by 

A(z' ■ = 52{-27rA^!J^,(5cosr 

+ sin2T + y - a/3Ci 2'sin^ r|. (33) 

Despite an early measurement of 1 nC /cm^ Q , the mag¬ 
nitude of the polarization change A(z' • below He 

extrapolated to zero field has recently been estimated as 
40 nC/cm^ pQ. 

The Hamiltonian in zero electric field can be simply 
written in terms of the induced polarizations as 


+ x(l) ^ -^1 

^z'v' 


sc 


pMS _ 


Al?) 






"z' 


>MS 


p“_ A,' 


>ANI 


2iibH 

N 


y]m ■ Si. 


(34) 


Introducing the field dependence of the DM interactions, 
we find = —y/2dDi/dEz' and = —dD 2 /dEz'. 
Similarly, Ciz' = V^dJi/dEz', C^z' = V^dJ 2 /dEz', 
and ^4 = AdKjdEz'. 

All z' components of the induced polarization ap¬ 
pear in H above. Because P^^ appears above Tn in the 
paramagnetic phase, each static magnetically-induced 
polarization along z' corresponds to a term in the Hamil¬ 
tonian. Due to the symmetry lowering associated with 
P^^, each bilinear spin term that appears in H also con¬ 
tributes to an induced polarization parallel to P^^. 

Taking c = d in Eq. ([TT|) , 


/pind^ _ 1 W 

^ ' AT 5E 


(35) 


has no components perpendicular to z'. Components 
of the operator P^^*^ perpendicular to z' would then 
contribute only to the transition matrix elements (n ^ 
0 |pind|o) other words, H includes all induced po¬ 
larizations with static contributions (0|P^^^|0) but not 
induced polarizations with only dynamical contributions 
in ^ OlP^^^IO). For example, P^^ does not appear in H 
because (P^^^ = 0. 

We used Eq. (|34|) to check our numerical results for the 
matrix elements (n|z' • P|0). Since {n\l-L\^) = the 

appropriate sum of polarization matrix elements with the 
field-dependent term — A/'id(n|m-M|0) must vanish when 
n 7 ^ 0. We verified that this condition is indeed satisfied 
for all excited states and magnetic fields. 


V. THZ ABSORPTION 


The absorption of THz light is given by (a(cd) = 
(2co’/c) ImA(co’) where^i^ 

~ ± vT(w) (36) 


is the complex refractive index for a linearly polarized 
beam, and y"^® are the dielectric, magnetic, 

and magnetoelectric susceptibility tensors describing the 
dynamical response of the spin systemp^~— and is 
the background dielectric constant tensor associated with 
charge excitations at higher energies. Subscripts i and 
j refer to the electric and magnetic polarization direc¬ 
tions, respectively. The second term, which depends on 
the light propagation direction and produces NDD, is 
separated from the mean absorption by writing N{uj) = 
N{oj)±xTi‘^iio). 

Summing over the SW modes n at the cycloidal order¬ 
ing wavevector Q, Aa{uj) = (dcj/c) Imy^®(co’) is given 
by 


AQ;(a;) = y] (5(w - w„), 

n 

(37) 


(38) 

po„ = (0|P-d-e/V|n), 

(39) 

HOn = (0|M ■ h/i2B\n), 

(40) 


where V = is the volume per Fe site, P^^^/V is given 
in units of nC/cm^ and 

_ 47r/iB nC _ 0.1388 
h cm^ cm 

The THz electric and magnetic fields are polarized in the 
e and h directions, respectively. 

After expanding N{uj) for small susceptibilities, we find 
that a{uj) = {2uj/c) ImN{uj) is given by 

a(,u)) ='^Bn5{u) - Un), (42) 


= Aw„|Fi|/9o„P + y2|M0np}, (43) 

where 

,, ttV nC^ 6.975 x 10“^ 

Yi = - ^—7 = -^-, (44) 

hc€o \/e“ cm^ \/e“ cm 

^ ^ 1-727v^ 

hcV cm 

Notice that X = Ay/Y\^. 

The dielectric constant depends on the polarization 
e of light. Based on a fit to the interference fringes, 
= 27.54 and 51.55 for e = [1,1,0] and e = [1,-1,0], 
respectively. 

For each orientation of the static magnetic field and 
light polarization, the integrated weight of every spec¬ 
troscopic peak at ujn is compared with the measured val¬ 
ues. This eliminates estimates of the individual peak 
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TABLE II: Fitting parameters (nC/cm^) from Aa 


X 


( 1 ) 


z'y' 




A 


( 2 ) 


iVn, 


Xmin 


1.543 


fit 1 -82.0 -50.3 +35.2 +13.6 4 

error ±3.1 ±8.0 ±1.9 ±3.0 

fit 2 -78.7 -39.4 +33.7 +13.9 3 1.536 

error ±3.5 ±1.7 ±1.9 ±3.0 


widths. Because the polarization and magnetization ma¬ 
trix elements are generally complex with an arbitrary 
overall phase that differs for each mode n, we can choose 
{0\Pz'\n) to be real. Other magnetization and polariza¬ 
tion matrix elements for mode n are then either purely 
real or imaginary. Under reversal of the field orien¬ 
tation, our numerical results indicate that (0|M|n) ^ 
— (0|M|n)* and (n|P|0) ^ (n|P|0)*. It follows that the 
NDD vanishes for zero field. Our numerical results also 
indicate that the NDD should vanish^ for field directions 
[0,0,1] and [1,1,0]. 

Unfortunately, fitting results for the mean absorption 
a{uj) were markedly inferior to results for Aa{uj). This 
may be caused by uncertainty about the dielectric con¬ 
stants e^, which does not enter Aa{uj). Moreover, the 
measured difference between the absorption in positive 
and negative fields is much less prone to systematic ex¬ 
perimental error than a{uj). 

Experimental results for the NDD with field along m = 
[1, —1, 0] are plotted in Figs.5(a) and (b) for e = [1, —1, 0] 
and [1,1,0], respectively. For some modes, the NDD is 
strong enough that a{uj) is small for light traveling in 
one direction but large for light traveling in the other 
direction^. In particular, for ~ 15.5 cm“^ and e = 
[1,1,0], An = 0.67 cm“^ for an 8 T field along m = 
[1,—1,0] and light propagating along k = [0,0,1] while 
An =4.12 cm“^ when either m or k is reversed. 

Fits to the NDD are based on the plotted 2, 4, 6, 8, 
10, and 12 T data sets. For each field value with two 
light polarizations, we evaluate the integrated weights for 
the 8 modes and between 

roughly 12 and 35 cm“^. Hence, there are A^dat = 96 
data points for A(a(co’). Because the and 

pairs are nearly degenerate in small fields, the average 
predicted NDD of those pairs are compared with the mea¬ 
sured values of Aa{uj) in a 2 T field. 

Remarkably, the NDD for m = [1, —1, 0] is dominated 
by the two sets of SC polarizations Pf^ and Pf^ associ¬ 
ated with the DM interactions Di and D2, respectively. 
Neglecting MS and ANI leaves four polarization param¬ 
eters: A^^^ = X^y,\ and . The diagonal 

parameters A^^^, = —A^!^, = c — d (domain 2) are ob¬ 
tained from X^y,\, = — 2\/2c and A^!^, = —y/2d using the 
relations in Appendix B. 

Above about 12 T, agreement between the theoretical 
and experimental values for Aa begins to deteriorate. 
This failure may be caused by the avoided mode cross- 



FIG. 5: (Color online) The measured (a,b) and predicted (c,d) 
NDD for 2 to 12 T fields along m = [1, —1,0] and for e = 
[1, —1, 0] (a and c) or [1,1, 0] (b and d). Predictions are based 
on fit 2. Dashed curves in (c) and (d) are the predicted mode 
frequencies. 


ings between 12 and 16 T, which mix 


and 


^2^^ and are imperfectly captured by our model. Since 
each of those modes exhibits pronounced NDD, the fitted 
polarization parameters are very sensitive to the precise 
behavior of the modes at the avoided crossings. 

Based on the typical noise level for the absorption, we 
set the experimental uncertainty for the integrated NDD 
of each peak to a = 1 cm“^. The error bars for each 
polarization parameter are then obtained from the con- 
dition that x^/xLin increases l/(A^dat - iVpar - !)• 

The results for fit 1 in Table II with Apar = 4 free 
parameters indicate that for domain 2, A^!^, = — A^!^, ~ 
—6.6 ± 6.7 nC/cm^. So our results imply the absence of 
the diagonal terms in A^^^ and support the simplified form 
of Eq.([T2]) for the first SC-induced polarization. With a 

slightly smaller Xmin^ ^ with Apar = 3 takes A^!^, = 

Xi' = 0 = 2\^X'- 

Because the sample may point parallel or antipar¬ 
allel to [1,1,1], the overall sign of Aa{uj) and of the polar¬ 
ization parameters is ambiguous. According to Katsura 
et however, Xz' = — A^!^,/a/2 in Eq.([T2]) should be 

positive so that A^!^, < 0. This condition is used to fix 
the overall sign of the SC parameters in Table II. 

Results for fit 2 are plotted in Figs.5(c) and (d). Al¬ 
though it underestimates the NDD for ^2^^ when e = 
[1,1,0] (perhaps due to a small shift in the spectra a{uj) 
for positive and negative fields), this fit otherwise de- 
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V (cnr^) V (cnr^) 


FIG. 6: (Color online) The measured NDD for 2 to 12 T fields 
along m = [1,1,0] and for e = (a) [1,—1,0] or (b) [1,1,0]. 
Dashed curves are the predicted mode frequencies. 


scribes all the relevant features of the NDD with field 
orientation m = [1,—1,0]. Based on fit 2 parameters, 
the predicted SC polarizations (Pi^) ~ 7.0 z' nC/cm^ 
and (Pf^) ~ — 2 . 1 z' nC/cm^ point parallel and antipar¬ 
allel to P^^, respectively. Hence, the total SC-induced 
polarization points along P^^ with a magnitude of 4.9 
nC/cm^. From Eq.(|32]), the SC-induced polarization of 
the canted AF above He is given by —4.2z' nC/cm^ op¬ 
posite to P™. 

Since the MS-induced polarizations above and below 
He differ by — sin^ r ^ 6.9 x 10“^ this 

term can be ignored in Eq. (| 33|) . Using the LSDA-hU re¬ 
sult [40| ^4 ^ 110 nC/cm^, the change in (z'-P^^d) below 
He contains the ANI-induced contribution 5'^^ 4/8 ~ 86 
nC/cm^. Therefore, the total predicted change A(z' • 
^ 96 nC/cm^ in the induced polarization from 
above to below He is more than twice larger than the 
recent experimental estimates of 40 nC/cm^ . Alter¬ 

natively, fitting the experimental jump to Eq. (| 33|) gives 
^4 = 40 nC/cm^, 40% smaller than the LSDA+U predic¬ 
tion. 

A possible explanation for this discrepancy is that we 
have slightly overestimated D 2 and So- As mentioned 
above, taking D 2 < 0.079 meV or So < 0.185 would 
stabilize a canted AE phase with spins tilted out of the 
{x' plane. Because such a state would recoup some 
ANI energy, the predicted jump in the induced polar¬ 
ization at He would be reduced. The planar canted AE 
phase of Eq.(j9|) would then appear above some higher 
critical field > He- Due to the non-coplanar AE 
phase, the lower AE mode a in Eig.4 would decrease with 
field between He and i7', vanish at i7', and increase with 
field only above i7'. 

Using = -^/2dDildE,, and = -dD^ldE,,, 

the results of fit 2 for A^tt and }& can be used to 
evaluate the dependence of Di and D 2 on an electric 
field applied along z'. Raman measurementsi^ indicate 
that the spectroscopic modes exhibit significant depen¬ 
dence on an electric field of 75 kV/cm along [0,1,0]. 
Eor an electric field of 100 kV/cm along [1,1,1], we find 


ADi/Di = 6.0 X 10-^ and AD 2 /D 2 = -6.4 x 10“^ Al¬ 
though very small, the change in Di will slightly increase 
the size of S and reduce the period of the cycloid. The 
change in D 2 will slightly reduce the tilt angle r. 


VI. DISCUSSION 

Although the distorted cycloid of BiEeOs is produced 
by the competition between magnetic interactions, the 
SC polarization dominates the NDD of BiEeOs. The 
NDD of BiEeOs along m = [1,—1,0] is well described 
by our model. Due to the pronounced NDD for when 
e = [ 1 , 1 , 0 ], BiEeOs laay be used as an optical diode, 
transparent to light traveling in one direction but opaque 
for light traveling in the opposite direction. Despite the 
successes of this model, several issues must be addressed. 

Eor light propagating along k = [0,0,±1], symme¬ 
try arguments^ and our numerical results indicate that 
NDD should be absent for m = with stable do¬ 
main 1 if either e or h coincides with Even for 

m = [ 1 , 1 , 1 ], where all three domains are degenerate^, 
the NDD should vanish if domains 2 and 3 are equally 
populated. While NDD is not observed for m = [0, 0,1], 
the NDD for m = [1,1,0] is plotted in Eig. 6 . The most 
pronounced NDD is seen near the avoided mode crossing 
of and around 8 T. 

To estimate the relative sizes of the NDD in fields along 
[1,1,0] and [1, —1, 0], we calculate the net squared NDD, 

(«) 

n 

where the sum runs over all modes between 2 and 12 T 
and An was defined by Eq. (| 37|) . Since U = 3.50 and 9.45, 
respectively, the observed NDD is substantially weaker 
for [1,1,0] than for [1, —1,0]. Because the contributions 
from metastable domains 2 and 3 cancel each other, they 
can not explain the NDD observed for m = [1,1,0]. 
While a population imbalance between metastable do¬ 
mains 2 and 3 would produce very weak NDD, domain 
1 is expected to predominate above a few Tesla. Mis¬ 
alignment of the crystal could produce the observed NDD 
when m = [ 1 , 1 , 0 ] but the excellent agreement between 
the measured and predicted mode spectrum in Eig.4(b) 
suggests that the sample is aligned quite well. 

The NDD for m = [1,1,0] probably arises from an 
optical misalignmenfp^ with the polarization vectors e 
and h rotated about k = [0,0,1]. Eor m = [ 77 , 77 ,/^], 
rt{a) ^ U( 7 r/ 4 ) sin^( 2 Qf), where a is the angle between 
e and [ 1 , — 1 , 0 ]. As shown in Eig.7 for m = [ 1 , 1 , 0 ] and 
[0,0,1], U((a) peaks oX a = 7 r/ 4 , i.e. when e = [1,0,0] 
and h = [0,1,0] or e = [0,1,0] and h = [—1,0,0]. Eor 
m = [ 1 ,— 1 , 0 ], Vt{a) ^ U( 0 ) — (D( 0 )/ 2 ) sin^( 2 (a) is pre¬ 
dicted to drop to about D( 0)/2 at a = 7 r/ 4 . Note that the 
squared NDD for an individual mode does not obey these 
relations: they are obeyed only by the net squared NDD 
summed over all modes. Measuring the NDD while ro- 
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FIG. 7: (Color online) The predicted Q versus a for fields 
along [1,—1,0] (red), [1,1,0] (blue), and [0,0,1] (green). SC 
parameters are obtained from fit 2. Inset shows the rotation 
of e and h about [0, 0,1]. 


tating the THz polarization vectors about [0,0,1] would 
help resolve questions about the NDD when m = [1,1,0]. 

While the SC dominates the dynamical response of 
BiFeOs, MS dominates its static properties. As argued 
elsewhere^, —^/3S‘^Ciz'Z^ dominates the in¬ 

duced polarization and opposes below Tn, in agree¬ 
ment with recent measurements^. A unified model of 
ferroelectricity^i^ also concludes that SC and ANI make 
minor contributions to (P^^^) compared to MS. 

The distinction between static and dynamics proper¬ 
ties in BiFeOs is not surprising. Since spin fluctuations 
6Si are transverse to the almost collinear, cycloidal spin 
state (Si), we find that x (Sj) ^ 0 but (5Si • (Sj) ^ 0 
for nearby sites i and j. Because the ANI is BiFeOs 
is extremely weak, spin fluctuations more strongly af¬ 
fect the SC-induced polarization than the MS- and ANI- 
induced polarizations. By contrast, the almost collinear 
spin structure of BiFeOs efficiently produces a static po¬ 
larization through the MS and ANI but not through the 
SC since (Si)-(S_^-) 7^ 0 and (l/N) ^^{Siz'){Siz') 7^ 0 but 
(Si) X (S 7 «0. 

Tokunaga et al^ recently attributed the induced 
transverse polarization along y' to the first SC polariza¬ 
tion P?^ with matrix elements and Those 

authors found that |Af3L| ~ 104 nC/cm^ and lA^Ll « 73 
nC/cm^. By contrast, fit 1 indicates that |A^!^/| ^ 6 
nC/cm^ is very small. The result |A^!^/| from Ref. na is 

reasonably close to the result |A^!^/| ^ 50 ± 8 nC/cm^ 
from fit 1. 

Considering only the first set of SC terms asso¬ 


ciated with earlier work^i identified as an 

electromagnon^i^ that can be excited by a THz elec¬ 
tric field when H = 0. When both sets of SC terms are 
considered, and also become electrically 

active at zero field. Of the four modes observed in zero 
field, only at 20.4 cm“^ is not electrically active. 
Using the SC parameters in Table II, couples most 
strongly of all modes to a THz electric field for domains 
2 and 3. This mode also exhibits the strongest NDD for 
nonzero field. 

To summarize, the SC polarization matrix elements 
dominate the NDD in BiFeOs. But work remains to un¬ 
derstand the origin of the MS-induced and perpendicular 
polarizations in this important material. Our explana¬ 
tions for the jump in the induced polarization at He and 
for the observed weak NDD when m = [1,1,0] need to 
be confirmed. Nevertheless, we believe that the present 
work on BiFeOs provides a compelling example of how a 
quantitative microscopic theory of magnetoelectric cou¬ 
plings follows from an analysis of the observed dynam¬ 
ical magnetoelectric response based on an effective spin 
model supplemented by first-principles calculations. 
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and UN acknowledge support by the Estonian Ministry 
of Education and Research Grant IUT23-03 and by the 
Estonian Science Eoundation Grant ETE8703. 


Appendix A: Experimental Details 

A single ferroelectric domain BiEeOs sample with face 
area 7 mm^ (0,0,1) and thickness d = 0.37mm along 
[0,0,1] was grown at Rutgers University. Voigt measure¬ 
ments (k T m) up to 17 T were performed in Tallinn. 
Earaday (k || m) measurements up to 12 T were per¬ 
formed in Tallinn and up to 31 T in Nijmegen, as re¬ 
ported earlier—. 

The Tallinn laboratory uses a Martin-Puplett type in¬ 
terferometer with a Si bolometer operating at 0.3 K and 
a mercury arc light source. Light pipes direct light to 
the sample in a He exchange gas-filled sample chamber 
placed into the cold 52 mm bore of a vertical-field super¬ 
conducting 17 T solenoid. In the Voigt configuration, 
mirrors before and after the sample change the light di¬ 
rection perpendicular to m. A rotatable wire grid on 
the dielectric substrate polarizer is placed before the first 
mirror. The sample can be rotated about the axis paral¬ 
lel to the direction of light propagation. A set of low pass 
filters with different cut-off frequencies is situated on the 
filter wheel in liquid He between the sample chamber and 
the bolometer chamber below the solenoid. 





Applying a 17 T field at 4 K for tens of minutes pop¬ 
ulates magnetic domain 1 when m = [ 1 , 1 , 0 ] or domains 
2 and 3 when m = [1,—1,0] M- Spectra were then 
measured in different ±H fields for about 15 minutes per 
field. No change in the magnetic domain populations was 
observed when a —17 T field was applied after a +17 T 
field. 

The zero-field absorption spectrum was subtracted 
from the spectra measured in field, thereby canceling out 
diffraction and interference effects caused by the sample. 
The differential absorption coefficient is a{H) — (a(0) = 
— ln(/i:///o)/d, where Iq and Ih are light intensity spec¬ 
tra in zero and H field and d is the sample thickness. 
Negative peaks in the differential absorption spectra for 
all field values were used to calculate the zero-field spec¬ 
trum. To generate the field-dependent spectra, the calcu¬ 
lated zero-field spectrum was aNDDed to the differential 
spectra. The NDD spectra Aa = (a(H) — (a(—H) do not 
depend on the zero-field spectra. 


Appendix B: SC-induced polarizations 


The Di term in the Hamiltonian can be written 


1 /SC ^ p(«) . (g. X sy. 

(Bl) 



In the absence of an electric field. 


fP) = ^( 0 , 1 ,- 1 ), 

(B2) 

F(J') = ^(-1,0,1), 

(B3) 

F("^ = ^(r-i,o) 

(B4) 


along X, y, and z, respectively. The polarization associ¬ 
ated with Di is then 

sc 1 

N dEa 

= V E nS(SixS,)^, (B5) 

where = —dF^^/dE^. LSDA+[/ calculations^ re¬ 
veal that 

/O -d d\ 

= 0 -c -c , (B 6 ) 

\0 c c ) 


/ c 0 c \ 

= d 0 -d , (B7) 

\-c 0 -c ) 


n(^) 
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—c —c 0 \ 
c c 0 j . 

-d d 0 / 

(B 8 ) 


Consequently, is not parallel to 

In the lab reference frame {x^y^z}, regrouping terms 
for domain 2 yields with A^^^ = 

n(^) _ or 


( c c-d d \ 

—c —2c —c I . (B9) 

d c — d c J 


We transform this matrix into the cycloidal reference 
frame {x'^y'^z'} using the unitary matrix U for domain 
2 : 

/ 1/V2 0 -1/V2\ 

u= -i/Ve +73 -i/ve (BIO) 

V 1/V3 1/+ 1/V3 / 

SO that 

yi) = uA^^'>u-^ 

/c-d 0 0 \ 

= 0 d-c - 2+0 . (Bll) 

\ 0 -V2d 0 / 

The diagonal terms change sign for domains 1 and 3. 
The D 2 term in the Hamiltonian can be written 
= D 2 N z' • T 2 . Thus, the SC-induced polariza¬ 
tion associated with D 2 can be generally written = 
S /3 ^+3 ^ 2 / 3 - In the lab reference frame, A^^^ is given 

by^ 


A(2) 



(B12) 


Transforming into the cycloidal reference frame. 


(e-f 0 0 \ 

A^ 2 ) _ 0 e-f 0 (B13) 

VO 0 e + 2 / y 


for all three domains. So is diagonal with compo¬ 
nents A^^^ = = e — f and A^?^ = e + 2 /. 


Appendix C: MS-induced polarizations 


The MS-induced polarizations are tV"^ Wi^ 

and P 2 + = According to LSDA+f/ 

calculations^, are given in the lab reference frame 

by 

( g h h\ 

yi) = \ h g h\, (Cl) 

\h h g 
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p(2) 


0 -j j \ 

3 0 -j \ , 

-j 3 0 J 


(C 2 ) 


where j = g — h. Transforming into the cycloidal refer¬ 
ence frame, 


/j 0 0 \ 

= 0 J 0 , (C3) 

\ 0 0 g + 2h J 


( 0 -x/3i 0\ 

^£(2)^-1= 0 0 (C4) 

V 0 0 0 / 

for all three domains. It follows that Ciy' = j, Ciz' = 
g 2h, and C 2 = Therefore, Ciy' = C 2 /\/ 3 . 


Appendix D: ANI-induced polarizations 

The perpendicular ANI-induced polarization is 

associated with the dependence of the polarization direc¬ 
tion n on an electric field. The ANI energy is given by 
yANi _ —K^-{Si • n^)^, where 


Ui = [sm0izcos(l)iz,sm0izsm(l)iz,cos0iz] (Dl) 


dOizIdEz where Oix = cos • x) and Oiy = cos • 
y). Similar identities hold for (j)io,. It follows that 


pANI 


1 ^ 1 v- f 

N dE^ ~ ~N y'y dOi^ 'Ml 


I I 5£ANI qj^ 

^ d4>i^ dE^] ~ N dK 


1 dyANi _ 1 f dOiy 

N dEy “ “V dOiy ^ 


ayANi I 1 QYAm QJ^ 

d(t>iy dEy I N OK dEy' ^ ^ 

The first terms in Eqs. (ID2IID4I) then produce the = 
-Ay/QKdei^ /dEz polarization perpendicular to z'. Be¬ 
cause d(l)ia/dEa is modulated by (— 1 )’^% the second 
terms produce the ^3 = (2Kd(\)izldEz polar¬ 
ization perpendicular to z'. The final terms produce the 
^4 = AyJ^dK/dEz polarization along z'. 

In the lab reference frame, the perpendicular polariza¬ 
tions produced by the dependence of the polar and az¬ 
imuthal angles Oi and <pi on the electric field E are given 
respectively by 


EiSia - Sip)SY 


(D5) 


is the local single-ion ANI axis and Oiz = cos ^{ui • z). 
Consequently, 


pANI 


1 ^ 1 f QQ.^ 

N dE^ ~ ~N dOi, ^ 


51/ANI I ^ QyANl QJ^ 

^ d4>iz dE, J “ N dK dE,^ 


(D2) 


(pANi(2))^ = ^ ^ - Si^)Sis. (D 6 ) 

i,gjS 

In the cycloidal reference frame, these polarizations are 
given by Eqs.(j27j) and ([29]). 


which must be evaluated in the limit -A z' . Due to the 
rhombohedral crystal structure, dOixjdEx = dOiy/dEy = 


^ J.R. Teague, R. Gerson, and W. J. James, Solid State Com- 
mun. 8, 1073 (1970). 

^ I. Sosnowska, T. Peterlin-Neumaier, and E. Steichele, J. 

Phys. C: Solid State Phys. 15, 4835 (1982). 

^ D. Lebeugle, D. Colson, A. Forget, and M. Viret, Appl. 

Phys. Lett. 91, 022907 (2007). 

^ A.M. Kadomtseva, A.K. Zvezdin, Yu.F. Popv, A.P. Py- 
atakov, and G.P. Vorob’ev, JTEP Lett. 79, 571 (2004). 

^ M. Tokunaga, M. Azuma, and Y. Shimakawa, J. Phys. Soc. 
Jpn. 79, 064713 (2010). 

® J. Park, S.-H. Lee, S. Lee, F. Gozzo, H. Kimura, Y. Noda, 
Y.J. Choi, V. Kiryukhin, S.-W. Cheong, Y. Jo, E.S. Choi, 
L. Balicas, G.S. Jeon, and J.-G. Park, J. Phys. Soc. Jpn. 
80, 114714 (2011). 


^ S. Lee, M.T. Fernandez-Diaz, H. Kimura, Y. Noda, D.T. 
Adroja, S. Lee, J. Park, V. Kiryukhin, S.-W. Cheong, M. 
Mostovoy, and J.-G. Park, Phys. Rev. B 88, 060103 (2013). 
® D. Lebeugle, D. Colson, A. Forget, M. Viret, A.M. Bataille, 
and A. Gukasov, Phys. Rev. Lett. 100, 227602 (2008). 

^ M. Ramazanoglu, W. Ratcliff II, Y.J. Choi, S. Lee, S.- 
W. Cheong, and V. Kiryukhin, Phys. Rev. B 83, 174434 
( 2011 ). 

J. Herrero-Albillos, G. Catalan, J.A. Rodriguez- 
Velamazan, M. Viret, D. Colson, and J.F. Scott, J. 
Phys.: Condens. Matter 22, 256001 (2010). 

1. Sosnowska and R. Przenioslo, Phys. Rev. B 84, 144404 

( 2011 ). 

S. Lee, W.M. Ratcliff II, S.-W. Cheong, and V. Kiryukhin, 

















14 


AppL Phys. Lett. 92, 192906 (2008); S. Lee. T. Choi, W. 
Ratcliff II, R. Erwin, S.-W. Cheong, and V. Kiryukhin, 
Phys. Rev. B 78, lOOlOl(R) (2008). 

J. Jeong, E.A. Goremychkin, T. Guidi, K. Nakajima, G.S. 
Jeon, S.-A. Kim, S. Eurukawa, Y.B. Kim, S. Lee, V. 
Kiryukhin, S.-W. Cheong, and J.-G. Park, Phys. Rev. Lett. 
108, 077202 (2012). 

M. Matsuda, R.S. Eishman, T. Hong, C.H. Lee, T. 
Ushiyama, Y. Yanagisawa, Y. Tomioka, and T. Ito, Phys. 
Rev. Lett. 109, 067205 (2012). 

Z. Xu, J. Wen, T. Berlijn, P.M. Gehring, C. Stock, M.B. 
Stone, W. Ku, G. Gu, S.M. Shapiro, R.J. Birgeneau, and 

G. Xu, Phys. Rev. B 86, 174419 (2012). 

M. Cazayous, Y. Gallais, A. Sacuto, R. deSousa, D. Lebeu- 
gle, and D. Colson, Phys. Rev. Lett. 101, 037601 (2008). 
P. Rovillain, R. de Sousa, Y. Gallais, A. Sacuto, M.A. 
Measson, D. Colson, A. Eorget, M. Bibes, A. Barthelemy. 
and M. Cazayous, Nat. Mat. 9, 975 (2010). 

D. Talbayev, S.A. Trugman, S. Lee, H.T. Yi, S.-W. 
Cheong, and A.J. Taylor, Phys. Rev. B 83, 094403 (2011). 
U. Nagel, R.S. Eishman, T. Katuwal, H. Engelkamp, D. 
Talbayev, H.T. Yi, S.-W. Cheong, and T. Room, Phys. 
Rev. Lett. 110, 257201 (2013). 

I. Sosnowska and A.K. Zvezdin, J. Magn. Magn. Mat. 140- 
144, 167 (1995). 

A.P. Pyatakov and A.K. Zvezdin, Eur. Phys. J. B 71, 419 
(2009). 

K. Ohoyama, S. Lee, S. Yoshii, Y. Narumi, T. Morioka, H. 
Nojiri, G.S. Jeon, S.-W. Cheong, and J.-G. Park, J. Phys. 
Soc. Jpn. 80, 125001 (2011). 

C. Ederer and N.A. Spaldin, Phys. Rev. B 71, 060401(R) 
(2005). 

M. Ramazanoglu, M. Laver, W. Ratcliff H, S.M. Wat¬ 
son, W.C. Chen, A. Jackson, K. Kothapalli, S. Lee, S.-W. 
Cheong, and V. Kiryukhin, Phys. Rev. Lett. 107, 207206 
( 2011 ). 

In previous work^^, the SW frequencies were scaled by 
\JS{S + 1). This artificial introduction of quantum fluc¬ 
tuations made it awkward to obtain the correct SW 
frequencies of the canted AE phase, where the upper 
mode frequency is proportional to 2fiBH rather than 
2jjjB^{S + 1)/SH. In the present work, the frequencies are 
scaled by S. Consequently, the Heisenberg interactions in¬ 
crease from Ji = —4.50 meV and J 2 = —0.20 meV to 
Ji = —5.32 meV and J 2 = —0.24 meV. This rescaling was 
also applied to the quoted values of Di and K in Ref. . 
Also keep in mind that we have redefined Di to be \/2 
larger than in earlier work. 

J. Jeong, M.D. Le, P. Bourges, S. Petit, S. Eurukawa, S.- 
A. Kim, S. Lee, S.-W. Cheong, and J.-G. Park, Phys. Rev. 
Lett. 113, 107202 (2014). 

R.S. Eishman, J.T. Haraldsen, N. Eurukawa, and S. Miya- 
hara, Phys. Rev. B 87, 134416 (2013). 

R. S. Eishman, Phys. Rev. B 87, 224419 (2013). 

J.J. Hopfield and D.G. Thomas, Phys. Rev. Lett. 4, 357 
(1960). 

D. Szaller, S. Bordacs, 1. Kezsmarki, Phys. Rev. B 87 
014421 (2013). 

1. Kezsmarki, N. Kida, H. Murakawa, S. Bordacs, Y. 
Onose, and Y. Tokura, Phys. Rev. Lett. 106, 057403 
( 2011 ). 

S. Bordacs, 1. Kezsmarki, D. Szaller, L. Demko, N. Kida, 

H. Murakawa, Y. Onose, R. Shimano, T. Room, U. Nagel, 
S. Miyahara, N. Eurukawa, and Y. Tokura, Nat. Phys. 8, 


734 (2012). 

S. Miyahara and N. Eurukawa, J. Phys. Soc. Japan 80, 
073708 (2011). 

1. Kezsmarki, D. Szaller, S. Bordacs, B. Kocsis, Y. Toku- 
naga, Y. Takuchi, H. Murukawa, Y. Tokura, H. En¬ 
gelkamp, T. Room, and U. Nagel, Nat. Comm. 5, 3023 
(2014). 

Y. Takahashi, Y. Yamasaki, and Y. Tokura, Phys. Rev. 
Lett. Ill, 037204 (2013). 

36 Takahashi, R. Shimano, Y. Kaneko, H. Murakawa, and 
Y. Tokura, Nat. Phys. 8, 121 (2012). 

I. Kezsmarki, U. Nagel, S. Bordacs, R.S. Eishman, J.H. 
Lee, H.T. Yi, S.-W. Cheong, and T. Room (unpublished). 
H. Katsura, N. Nagaosa, and A.V. Balatsky, Phys. Rev. 
Lett. 95, 057205 (2005). 

R.S. Eishman, Phys. Rev. B 88, 104419 (2013). 

J. -H. Lee and R.S. Eishman, (unpublished). 

R. deSousa and J.E. Moore, Phys. Rev. B 77, 012406 
(2008). 

When m = [0 , 0. 1], a field-independent feature was re¬ 
ported in Ref. [1^ at 16.3 cm“^. Because it did not fit 
into our mode scheme, this feature was difficult to under¬ 
stand. A re-examination of the THz data reveals that the 
total area of the nearly-degenerate modes at zero 

field was erroneously reported too high (6.9± 1.0 cm“^) by 
Ref. 0. The correct value is 3.5 ± 0.3 cm ^ and there is 
no magnetic field independent mode at 16.3 cm“^. Hence, 
it does not appear in Eig.4(a). 

It is questionable whether Ref. [3 actually observed the 
nearly-degenerate modes at zero field: the re¬ 

ported frequency of 18.4 cm“^ is substantially higher than 
the 16.2 cm“^ mode seen in Eig.4. In Earaday measure¬ 
ments, Ref. 0 also observed a feature in the absolute spec¬ 
trum around 18 cm“^ at zero field. Since this feature does 
not appear in a magnetic field, it may be an interference 
fringe rather than an actual mode of the cycloid. 

C. Jia, S. Onoda, N. Nagaosa, and J.H. Han, Phys. Rev. 
B 76, 144424 (2007) 

T. A. Kaplan and S.D. Mahanti, Phys. Rev. B 83, 174432 

( 2011 ). 

M. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006). 

LA. Sergienko and E. Dagotto, Phys. Rev. B 73, 094434 
(2006). 

R. deSousa, M. Allen, and M. Cazayous, Phys. Rev. Lett. 
110, 267202 (2013). 

A polar vector does not change sign under spatial inversion 
but an axial vector does. Under spatial inversion, 

m,. =i s,.s, 

,Rj=Ri + ax 

^ s 2: = 

R^ ,Rj =Ri —ax 

Since similar relations hold for Wi^y and Wi^z, Wi is a 
polar vector. On the other hand, 

X (-irsi-s,- 

Rj ,Rj=Rj+ax 

^ V ^ (-irsi-s, = -H/2,=c 

R^ ,Rj=Rj —ax 

because Ui = + 1 for and Hj = — ax. Therefore, 

W 2 is an axial vector but z' x W 2 is a polar vector. 




15 


S. Miyahara and N. Fumkawa, Phys. Rev. B 89, 195145 
(2014). 

Assuming the simple coplanar cycloid of Eqs. (IBISli . symme¬ 
try dictates (I. Kezsmarki, unpublished) that Aa(cL;) = 0 
when m J_ and either e or h lies along for do¬ 
main m. Those two conditions are met for domain 1 with 
Xi = [ 1 ,— 1 , 0 ] when m = [ 77 , 77 ,/c], k = [ 0 , 0 , ± 1 ], and ei¬ 
ther e = [lA^O] or [1,—1,0]. This argument implies that 
NDD should vanish for static magnetic fields along [0, 0,1] 
and [ 1 , 1 , 0 ] but appear for held along [ 1 ,— 1 , 0 ], when do¬ 
mains 2 and 3 are stable and degenerate^. Although this 
symmetry argument fails when the cycloid is distorted by 
a magnetic held, our numerical results suggest that AQ;(ct;) 
vanishes when the two conditions above are satished. 

See, for example, P.R. Bevington, Data Reduction and Er¬ 
ror Analysis for the Physical Sciences (McGraw Hill, New 
York, 1969), Section 11.5. The parameters in Table II for ht 
2 are very close to those in Ref. [13] , which did not include 
the 2 T data set due to its small NDD. The error bars in 


Table II are also close to those in Ref. [13 , which undertook 
a more sophisticated analysis based on the covariance. 
About 3% of the light intensity has the wrong (perpen¬ 
dicular) polarization. This misalignment is more severe for 
the Voigt geometries (m = [1,—1,0] and [1,1?0]) than for 
the Faraday geometry (m = [0, 0,1]). 

J.-H. Lee, I. Kezsmarki, and R.S. Fishman, (unpublished). 
H.J. Xiang, P.S. Wang, M.-H. Whangbo, and X.G. Gong, 
Phys. Rev. B 88 , 054404 (2013). 

X.Z. Lu, X. Wu, and H.J. Xiang, Phys. Rev. B 91, 
100405(R) (2015). 

M. Tokunaga, M. Akaki, T. Ito, S. Miyahara, A. Miyake, 
H. Kuwahara, and N. Fumkawa, Nat. Gomm. 6 , 5878 
(2015). 

S. Miyahara and N. Fumkawa, J. Phys. Soc. Japan 81, 
023712 (2012). 

H.B. Ghen and Y.-Q. Li, Applied Physics Letters 102 , 
252906 (2013). 


